Constraints on inelastic dark matter from XENON10 
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It has been suggested that dark matter particles which scatter inelastically from detector target 
nuclei could explain the apparent incompatibility of the DAMA modulation signal (interpreted as 
evidence for particle dark matter) with the null results from CDMS-II and XENON10. Among 
the predictions of inelastically interacting dark matter are a suppression of low-energy events, and 
a population of nuclear recoil events at higher nuclear recoil equivalent energies. This is in stark 
contrast to the well-known expectation of a falling exponential spectrum for the case of elastic 
interactions. We present a new analysis of XENON10 dark matter search data extending to E nr = 
75 keV nuclear recoil equivalent energy. Our results exclude a significant region of previously allowed 
parameter space in the model of inelastically interacting dark matter. In particular, it is found that 
dark matter particle masses m x > 150 GeV are disfavored. 



PACS numbers: 95.35. +d, 14.80. Ly, 29.40.Gx, 95.55.Vj 



I. INTRODUCTION 



Astrophysical evidence indicates that 23% of the mass 
of the universe is in the form of non-baryonic dark matter 
[TH5]. A well- motivated dark matter candidate particle 
is found in supersymmetric extensions to the Standard 
Model, in which the lightest supersymmetric particle 
(LSP) is stable. A cosmologically interesting relic density 
of Weakly Interacting Massive Particles (WIMPs) arises 
from rather general arguments 0] , which implies that the 
LSP could be dark matter. The open question of the ex- 
pected mass and cross-section of particle dark matter is 
being addressed by numerous direct and indirect detec- 
tion experiments E] . The nature of the dark matter 
particle, and its coupling to standard model particles has 
been the subject of much recent theoretical investigation 

EH2]. 
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A. The XENON10 detector 

The XENON10 detector is a liquid xenon time- 
projection chamber with an active target mass of 13.7 kg. 
It was designed to directly detect galactic dark matter 
particles which scatter off xenon nuclei. Typical veloc- 
ities of halo-bound dark matter particles are of order 
10~ 3 c. This leads to a a featureless exponential recoil 
energy spectra for spin-independent elastic scattering of 
dark matter particles on a xenon target. In the case of a 
100 GeV dark matter particle mass, the predicted elas- 
tic scattering recoil energy spectrum falls by an order 
of magnitude from — 30 keV nuclear recoil equivalent 
energy (keVr) (IHU). 

A particle interaction in liquid xenon creates both ex- 
cited and ionized xenon atoms [11], which react with 
the surrounding xenon atoms to form excimers. The ex- 
cimers relax on a scale of 10 -8 s with the release of scintil- 
lation photons. This prompt scintillation light is detected 
by 88 photo-multiplier tubes and is referred to as the 5*1 
signal. An external electric field {E c i = 0.73 kV/cm) 
across the liquid xenon target causes a large fraction of 
ionized electrons to be drifted away from an interaction 
site. The electrons are extracted from the liquid xenon 



2 



and accelerated through a few mm of xenon gas by a 
stronger electric field (~ 10 kV/cm), creating a secondary 
scintillation signal. This scintillation light is detected by 
the same photo- multiplier tubes, is proportional to the 
number of ionized electrons and is referred to as 5*2. 

The XENON10 detector discriminates between elec- 
tron recoil background and the expected nuclear recoil 
signal from the scattering of a dark matter particle, via 
the distinct ratio of proportional (52) to primary (51) 
scintillation for each type of interaction. The XENON10 
collaboration previously reported WIMP-nucleon exclu- 
sion limits for spin-independent [12] and spin-dependent 
[T3"] elastic scattering. 



B. Inelastic dark matter 

In stark contrast to elastic scattering jHEj], a model of 
inelastic interactions of dark matter (iDM) predicts [TU 
[T7] a suppression of low-energy nuclear recoils. The recoil 
energy spectrum peaks near 40 keVr, for a 100 GeV dark 
matter particle incident on a xenon target (and standard 
halo assumptions). As explained in [17] . this feature is 
a result of the mass difference S between the proposed 
ground and excited states of the dark matter particle. 
The value of S is unknown and is a free parameter in 
the model, subject to the physical constraint that 6 > 
170 kcV would kincmatically forbid any scattering from 
dark matter particles, given an expected particle mass 
m x = 100 GeV and a galactic escape velocity v esc = 
500 km s _1 . In the limit S — > 0, the iDM model reduces 
to the usual elastic scattering case. 

At present, the iDM model is comfortably consistent 
[T8"] with all reported null results (including the recent 
CDMS results [TS]) as well as the claimed detection from 
DAMA [H]. Interpretation of the XENON10 results in 
[18] was limited by the fact that in [12J . data were not 
published beyond 45 keVr. Since the focus in [T2] was 
on elastic interactions, that analysis was optimized for 
events with recoil energies in the range 4.5 — 26.9 keVr. 
This paper presents a new analysis of the XENON10 dark 
matter search data, extending to the energy range of in- 
terest for inelastic dark matter scattering. In doing so, 
we find new constraints on allowed parameter space in 
the iDM model. 



II. PREDICTED EVENT RATES 

We calculate predicted differential event rates as a 
function of nuclear recoil energy in XENON10 and 
DAMA following [18], as 

<W N T M N J^A*F\E nr ) r M dv. (1) 
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is A. We assume the standard local dark matter density 
p x = 0.3 GeV cm -3 , with dark matter particle mass m x 
and cross section (per nucleon) a n . Therefore, the re- 
duced mass n ne is for the nucleon — dark matter particle 
system. 

The nuclear form factor F(E nr ) accounts for a loss of 
coherence as momentum transfer to the nucleus increases. 
We use the Helm form factor parameterization F(E nr ) = 
3ji(qr n )/qr n ■ exp(— (qs) 2 /2). We take the effective nu- 
clear radius r n = \/r 2 — 5s 2 , with r = 1.2A 1 / 3 fm and 
the skin thickness s = 1 fm. The momentum transfer to 
the nucleus is just q = \/2M]yE nr . It was shown in [TU] 
that F 2 using this parameterization differs by 8% (13%) 
at 30 keVr (60 keVr) compared with other reasonable pa- 
rameterizations. As noted in [IS], this can lead to differ- 
ences of up to 25% in the predicted rate. For consistency, 
we therefore use this same form factor to calculate pa- 
rameter space consistent with DAMA and excluded by 
XENON10. 

We assume the dark matter velocity distribution to be 
Maxwellian, and perform the integration over f(v)/v fol- 
lowing [20 . As in [20] , we take the rotational speed of the 
local standard of rest to be vq = 220 km s _1 . We used 
a weighted average of the velocity distribution over the 
period October 2006 — » February 2007, corresponding to 
when the data were acquired. The lower limit of the in- 
tegral in Eq. [l]is x/l/ (2m N E nr )(m N E nr /n + S), 
which reduces to the usual v m i n (as in [10], for example) 
in the case S 0. The upper limit is set by the galactic 
escape velocity, which we take to be v esc = 500 km s _1 
as in |18j . Uncertainties in exclusion limits arising from 
assumptions about v esc are discussed further in Sec. |Vj 



III. EVENT ACCEPTANCE 

Figure [l] shows the results of a re-analysis of 58.6 live 
days of dark matter search data already described in [12l 
ITS"] , with the fiducial target mass unchanged at 5.4 kg. 
While the previous analyses were performed blind, the 
present work introduces two new software cuts (discussed 
in Sec. Ill A and Sec. Ill FJ I and so cannot be considered 



The number of target nuclei in the detector is Nt, the 
mass of the target nucleus is Mn and its atomic number 



as blinded. The hardware trigger efficiency was verified 
to be greater than 99% for 4 electrons in the 52 signal, 
and a software threshold of 12 electrons was imposed in 
the analysis. 

The recoil energy range is extended to the region rel- 
evant for inelastic dark matter scattering, and is quoted 
in keVr. The energy scale is calculated from the primary 
scintillation signal as E nr = Sl/(L y £ e //)-(5 e /5„), with 
L y = 3.0±0.14 photo-electrons/keVee the measured light 
yield for 122 keV photons. The scintillation quenching of 
electron and nuclear recoils due to the electric field E^ 
are 5 e = 0.54 ± 0.02 21j and 5„ = 0.95 ± 0.05 [33]. 
C e f f is the relative scintillation efficiency of liquid xenon 
for nuclear recoils. Due to the varied results from ex- 
perimental determinations of C e f f |21H25j , we show our 
results in Fig. [Tjassuming a constant C e ff — 0.19. From 
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this, we find 1.00 ±0.15 51 photo-electron/keVr. Prior 
to calculating exclusion limits, appropriate values of keVr 
for each remaining event were obtained based on recent 
measurements of C e f / 




30 40 
keVr 



FIG. 1: (color online) Events remaining after re-analysis of 
58.6 live days of dark matter search data. The dark matter 
acceptance box is bounded by the stair-step lines indicating 
the centroid and — 3cr bounds obtained from fits to the neu- 
tron calibration data. The fitted electron recoil centroid is 
shown as a solid blue curve, and the — 3<r contour is shown 
dashed. The 13 events remaining in the acceptance box are 
circled. An 52 software threshold of 12 electrons was imposed 
in the analysis. 

Background rejection is obtained for each event from 
the ratio of proportional to primary scintillation, written 
as log 10 (52/51). An acceptance box for nuclear recoils 
was defined from the nuclear recoil calibration data, fol- 
lowing the procedure described in |12) . This is indicated 
by the stair-step lines in Fig. [I] and the acceptance 
A nr is listed in Table |T] Additional background rejec- 
tion is accomplished with software cuts on 51 photo- 
multiplier tube optical distributions (described in Sec. 
Ill A I, and o n the 51 decay time distribution (described 
in Sec. 



IIIB) 



The 13 events from the dark matter search data that 
remained in the acceptance box after all cuts are indi- 
cated in Fig. [T] by red circles. No additional events re- 
mained in or below the acceptance box in the energy 
range 75 — 250 keVr, and in the interest of clarity this is 
not shown in Fig. [T] The acceptance remains A nr ~ 0.50 
in this higher energy range. The 5 events below the ac- 
ceptance box in Fig. [l]are consistent with the false single 
scatter pathology described in Sec. |IV| 



A. Poisson likelihood 

As explained in [12], a residual population of back- 
ground events were observed in the electron recoil cali- 
bration data. These events tend to populate the nuclear 



TABLE I: Software cut acceptance for nuclear recoils e p and 
e/ (discussed in Sec. Ill A and Sec. IIIB l, and the nuclear 
recoil band acceptance A nr , as a function of nuclear recoil en- 
ergy. The expected number of events N sta t in the acceptance 
box is determined from the number of detected events N e „t s 
and the Gaussian width of the electron recoil log 10 (52/51) 
distribution. 



E nr (keVr) e p ej A r . 



N s t a t 



02-05 
05-10 
10-15 
15-20 
20-25 
25-30 
30-35 
35-40 
40-45 
45-50 
50-55 
55-60 
60-65 
65-70 
70-75 



1.00 0.91 
0.92 0.90 
0.83 0.91 
0.67 0.89 
0.62 0.91 
0.59 0.92 
0.61 0.90 
0.64 0.92 
0.65 0.92 
0.56 0.92 
0.64 0.92 
0.63 0.90 
0.65 0.89 
0.63 0.89 
0.67 0.94 
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recoil signal region. The origin of these events is dis- 
cussed in Sec. |IV| They were targeted with high efficiency 
by calculating the Poisson likelihood of the 51 hit pat- 
tern of each event. For each event, the probability at that 
the i th photo-multiplier tube registered a photo-electron 
(resulting from the 51 primary scintillation signal) is a 
function of the event vertex, so on = Ui(x,y, z). As de- 
scribed in , the even t (x,y) coordinates are deter- 
mined by the pattern of proportional scintillation on the 
top photo-multiplier tube array, while the z coordinate is 
measured from the time delay between the primary and 
proportional scintillation. 



The distribution of cti for each photo-multiplier tube 
was measured from calibration data, obtained after in- 
troducing neutron-activated xenon into the XENON10 
detector. As described in [5S], this produced an internal, 
homogeneous source of 164 keV gamma rays from the de- 
excitation of 131m Xe. For each event, we then calculated 
the Poisson probability pi of obtaining the observed hit- 
pattern, given the expectation on. Two cut parameters 
were defined as Vb,t = log 10 (E a i/Pi)- In the case of 
'Pb CPt), the sum ran over only the bottom 41 (top 47) 
photo-multiplier tubes. Linear cut bands were optimized 
for Vb and Vt in order to remove leakage events in the 
electron recoil calibration data. This is discussed further 
in Sec. 



IV The acceptance for nuclear recoils e p for the 



combined Vb t software cuts is shown in Table |T| 
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B. Primary scintillation pulse shape 

It is well known that the primary scintillation light 
from nuclear recoils and electron recoils in liquid xenon 
exhibit distinct decay times corresponding to the pref- 
erentially excited singlet and triplet states of the Xe? 
dimer [27 . The resulting pulse shape discrimination 
is significantly less powerful than discrimination based 
on log 10 (52/51) (see [25], for example). This is due to 
the modest separation between the lifetime of the triplet 
(27 ns) and singlet (< 4 ns) states [29]. At small recoil 
energies, Poisson fluctuations in the number and arrival 
times of photo-electrons increase the width of the distri- 
bution, as shown in Fig. [2j It is also worth mention- 
ing that the 105 MHz clock speed of the ADC used in 
XENON10 is not optimal for characterizing xenon scin- 
tillation pulse shapes. 
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FIG. 2: (color online) (top) The distribution of f p for electron 
recoil calibration data. The black dashed line indicates the 
discrimination cut bound on this parameter, (bottom) The 
distribution of f p for nuclear recoil calibration data. The 
acceptance e/ of the cut is listed in Table [I] 

Still, a modest benefit can be obtained by parameter- 
izing the primary scintillation 51 = 51 (i) in terms of its 
prompt fraction (following [30]), 



) t t 'Sl{t)dt ' 



(2) 



The point at which a primary scintillation pulse rises 
to 10% of its maximum was defined as t , U = t — 5, 
if — to + 30 and the prompt window t w — 4. Times are 
quoted in units of 1/105 (is ADC samples, so t w w 38 ns. 
As discussed in [3D], the optimal value of t w depends on 
the electronics and the size of the detector. The discrim- 
ination was only slightly decreased if we chose t w = 3, 
while t w < 3 or t w > 4 showed significantly weaker per- 
formance. The prompt fraction f p obtained for neutron 
and electron recoil calibration data are shown in Fig. [2] 
along with curves indicating the centroid of each distri- 
bution. A cut which maintains e/ — 0.9 acceptance for 



nuclear recoils is indicated by the black dashed curve. 
The acceptance as a function of energy is shown in Table 
III 

The electron recoil rejection efficiency obtained from 
this cut rises monotonically from R er = 0.15 at E nr 
10 keVr to R er = 0.5 at E nr = 70 keVr. This level of re- 
jection is significantly weaker than the rejection already 
obtained from discrimination based on log 10 (52/51). 
However, it is an independent channel for characterizing 
the electromagnetic background, and is especially useful 
at higher energies. 



IV. ELECTROMAGNETIC BACKGROUNDS 

The signal acceptance box in Fig. [T] is subject to con- 
tamination ("leakage") from the electromagnetic back- 
ground. This background divides naturally into two com- 
ponents: statistical leakage, and (following the nomencla- 
ture of [15]) anomalous leakage. Predictions for the for- 
mer were obtained from Gaussian fits to the log 10 (52/51) 
distributions observed in calibration data, and are shown 
in Table|l](N stat ). Predictions for the latter, which will be 
shown to arise from mis-identified multiple scatter events, 
are not tabulated. This is due to the small number of 
such events in the calibration data sample. Instead, we 
calculate the predicted number of events for the energy 
range E nr < 75 keVr. 

Multiple scatter events in the active 13.7 kg region of 
the XENON10 detector are easy to identify, because each 
scatter vertex creates a separate proportional scintilla- 
tion (52) pulse. However, there are also 8.7 kg of xenon 
entirely outside the active region, and 1.2 kg of xenon 
in the region between the bottom photo-multiplier tube 
array and the cathode grid. A rendering of the detec- 
tor which shows these regions can be found in [24]. The 
1.2 kg of xenon below the cathode grid contribute 51 
but not 52, due to the reversed direction of the electric 
field Ed in this region. A multiple scatter event with 
one vertex in this region and one vertex in the central 
5.4 kg target would always — in the absence of the Vbt 
(or similar) characterization — be tagged as a valid single 
scatter. 

The 8.7 kg of xenon in the outer region were expected 
to be entirely passive, i.e. neither 51 nor 52 would be ob- 
served from energy depositions in this region. However, a 
ray-tracing Monte Carlo simulation [3T] showed that the 
light collection efficiency in this region ranged as high as 
15% (in a small region near the lower photo-multiplier 
tube array), with a mean of about 1%. 

Both the reverse field region and the outer region al- 
low the possibility that a multiple scatter background 
event is tagged as a single scatter event (a "false single 
scatter"). Some fraction of the 51 from the additional 
scatter may be collected, while none of the 52 will be col- 
lected. This implies that many false single scatters will 
have anomalously low values of log 10 (52/51), compared 
with genuine single scatter events. The lower bound in 
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log 10 (52/ 51) for false single scatter events is set by the 
52 software threshold, so these events may fall below the 
nuclear recoil acceptance box as well as in it. This is 
observed in electron recoil calibration data. 
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FIG. 3: (color online) The fraction of leakage events remain- 
ing among single gamma ray scatters in the electron recoil 
calibration data set (black stars). Also shown are the Monte 
Carlo prediction for false single scatters, with an additional 
scatter below the cathode grid (red triangles), or in the outer 
8.7 kg of xenon (blue circles). The sum of the two Monte 
Carlo components is indicated by pink squares. Note that 
the x axis is electron (not nuclear) recoil equivalent energy; 
20 keVee = 76 keVr. 

Figure [3] (black stars) shows the fraction of leakage 
events remaining among single scatters in an electron re- 
coil calibration data set, obtained using a 137 Cs source. 
This fraction was calculated prior to applying cuts based 
on Vb.t or f p , and only for events tagged as single scat- 
ters in the 5.4 kg target. An event was considered to be a 
leakage event if its log 1Q (S2/Sl) value fell below the — 3<r 
contour for electron recoils, which is shown as a dashed 
line in Fig. [I] To obtain the correct electron equivalent 
energy scale, the data were corrected for the measured 
scintillation yield as a function of energy [i32 . 

A Monte Carlo simulation of the XENON10 detector 
response to this electron recoil calibration was used to 
obtain a prediction for the false single scatter fraction. 
The result is shown in Fig. [3j divided into two com- 
ponents: false singles with an additional scatter below 
the cathode grid (red triangles), and false singles with 
an additional scatter in the outer 8.7 kg of xenon (blue 
circles). The sum of the two components is also shown 
(pink squares) . The simulation was corrected to account 
for the predicted light collection efficiency in each region 
[31]; for the scintillation quenching [21], based on the 
electric field strength in each region [3T]; and for the ob- 
served energy resolution of the detector. 

For each simulated false single scatter event, a value of 
log 10 (52/51) was calculated. This was done by assuming 
that each individual scatter vertex had a typical electron 
recoil log 10 (52/51) value, given by the centroid of the 
distribution (Fig. [I] solid blue line). The simulation 



prediction in Fig. [3] only includes false single scatters 
whose log 10 (52/51) was low enough that the event would 
be characterized as leakage. This criteria was met by 26% 
of events with an additional scatter below the cathode 
grid, and by 5% of events with an additional scatter in 
the outer 8.7 kg of passive xenon. 

It is clear from Fig. [3] that the observed leakage frac- 
tion and the predicted false single scatter fraction are 
very similar in the calibration data. Making this com- 
parison for the dark matter search data is significantly 
more involved (and is a work in progress |33j). since it 
requires detailed modeling of the location and isotopes 
which contribute to the electromagnetic background. 

From the electron recoil calibration data, the rejec- 
tion efficiency for false single scatter events in the nu- 
clear recoil acceptance box was calculated after all cuts 



(as described in Sec. Ill A and Sec. IIIB) to be R; ea fe = 
0.851 ± 0.056. If this same level of rejection applies to 
background events in the dark matter search data, there 
should remain 12.0 ±4.5 false single scatter events in the 
acceptance box in Fig. [TJ 

If we account for statistical leakage (N stat from Table 
|l| , the 13 remaining events in Fig. [ljwould be reduced to 
7.5 events. The rejection efficiency for the electron recoil 
calibration data in this case is Ki ea k = 0.908 ± 0.052. If 
this same level of rejection applies to background events 
in the dark matter search data, we predict 6.5±4.0 events 
remaining in the acceptance box in Fig. [T] In either 
case, the number of remaining potential signal events are 
consistent with expected background. 



V. DARK MATTER EXCLUSION LIMITS 

We calculate 90% CL exclusion limits on allowed re- 
gions of inelastic dark matter parameter space separately 
for a n — m x , a n — S and 5 — m x , treating all 13 remaining 
events in the acceptance box as potential dark matter 
signal. The results are shown in Fig. [4j for several rep- 
resentative cases. We use the pMax method from |34j . 
which explicitly provides for the presence of background 
events. This allows us to place conservative exclusion 
limits, without background subtraction (as in [35], for 
example). We verified that our results agree with [18] 
if we only consider events remaining in the acceptance 
box in [12] . take a constant e • A nr — 0.3, and a constant 
C e ff — 0.19. Regions of parameter space allowed by the 
DAMA modulation data [16] were calculated following 
the procedure described in [TS], and are shown as filled 
magenta (90% C.L.) and cyan (99% C.L.) regions. 

Our results strongly disfavor dark matter particle 
masses m x > 150 GeV, if we use the C e ff measurement 
from [2U [35] (Fig. [4j solid black curves). For lighter 
masses, our 90% CL exclusion limits still find common 
parameter space with the 99% CL DAMA-allowed region. 
The behavior of the exclusion limits are similar if we use 
the recent £ e ff measurement from [33] (Fig. |4j solid gray 
curves), however, larger dark matter particle masses re- 
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main consistent with DAMA for 110 keV < 6 < 130 keV. 
In this case, we also find that values of S < 100 keV are 
more strongly disfavored. 

The inelastic dark matter model is sensitive to the 
value of the galactic escape velocity, which we have 
fixed at v esc — 500 km s . If instead we take v esc — 
600 km s _1 , dark matter particles with higher kinetic 
energy would remain bound in the halo. The predicted 
event rate in XENON10 at higher recoil energies would 
thus increase, with the result that our data would more 
tightly constrain the DAMA-allowed parameter space. 

In addition to the exclusion limits shown in Fig. [4j it 
is worth mentioning that our results can also be viewed 
as a limit on the modulated fraction of the dark matter 
signal. A primary benefit of this quick analysis is that 
it requires no assumptions about the dark matter halo 
velocity distribution. The sole assumption is that the 
DAMA modulation is due to nuclear recoils from iodine 
nuclei, quenched relative to electron recoils by qi^ nT ~ 
0.09. Following [TB], the expected number of signal events 
observed by XENON10 can be related to the modulation 
amplitude observed by DAMA in the energy range 2 — 
6 keVee (22 — 67 keVr) according to 

N x = 4kS m (l/f - 1), (3) 

where the modulated amplitude is S m = 0.0131 ± 
0.0016 counts/kg/keVee as measured by DAMA [IB], and 
/ is the fraction of that signal that is modulated. The 
prefactor k = e • A nr (5A kg) (58.6 days) = 92 accounts 
for the total exposure obtained by XENON10. Account- 
ing for C e ff, 12 of the 13 remaining events fall in the 
specified energy range, so taking the 90% C.L. Poisson 
upper bound gives Nx — 18. From this, we find the 
model-independent result that a modulation / > 0.21 is 
required for the two experiments to remain consistent. 
Considering a 2er downward fluctuation in S m still re- 
quires / > 0.17. Such large modulation arises naturally 
for inelastic dark matter because of the minimum kinetic 
energy S that a dark matter particle must have in order 
to scatter. 

For completeness we note that the analysis presented 
here does not result in a significant change to the WIMP- 
nucleon exclusion limit quoted in |12j . for the case of 
spin-independent elastic dark matter particle scattering. 
The "maximum gap" method of |34j was used in |12j . 
Although the number of mis-identified events has de- 
creased dramatically in the present work, the maximum 



energy interval containing no potential signal events re- 
mains similar at approximately 4.5— 18 keVr, taking (the 
more conservative) C e ff measurement from |23j . 



VI. SUMMARY 

A primary conclusion of this work is that the events 
which have previously been referred to as "anomalous" 
background have been confirmed to originate in the pas- 
sive layers of LXe surrounding the target, as explained 
in Sec. |IV[ We have shown that these background events 
can be targeted with high efficiency by the cuts described 
in Sec. |III A| A corollary conclusion is that the occur- 
rence of this class of events can be explicitly prevented 
by the design of the detector, rather than by software 
cuts. Such a situation is realized if there is no xenon out- 
side the active target region, or (more realistically) if the 
photomultipliers viewing the target volume are perfectly 
blind to scintillation produced in passive areas of liquid 
xenon, including below the cathode grid. We have shown 
that the XENON10 dark matter search data exclude pre- 
viously allowed regions of the DAMA-allowed parameter 
space, in the model of inelastic dark matter. Specifi- 
cally, dark matter particle masses m x > 150 GeV are 
disfavored. While there are events remaining in the dark 
matter acceptance box in Fig. [T] and while these events 
appear consistent with the expected spectral shape for in- 
elastic dark matter, we do not claim a detection. It was 
demonstrated in Sec. |IV| that the number of remaining 
events is reasonably consistent with the expected num- 
ber of leakage events, based on electron recoil calibration 
data. 
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